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I. 


INTRODUCTION 


A.  WIRELESS  NETWORKS 

Wireless  Local  Area  Networks  (WLANs)  have  become  more  widespread  since 
the  ratification  of  the  IEEE  802.11b  standard.  WLANs  are  widely  deployed  in  places, 
such  as  office  rooms,  school  campuses,  companies,  and  even  bookstores  and 
coffeehouses  [1].  Unlike  the  wired  networks,  802.1 1-based  WLANs  broadcast  radio 
frequency  (RF)  data  and  provide  numerous  benefits,  such  as 

•  ease  of  installation, 

•  cost  effective, 

•  flexibility, 

•  mobility,  and 

•  improved  system  reliability. 

Some  network  managers  have  security  fears  regarding  WLANs  and  prefer  wired- 
LANs.  With  a  WLAN,  the  data  is  broadcast  between  the  client  devices  and  the  access 
points,  which  are  also  called  WLAN  endpoints.  Any  authorized  or  unauthorized  WLAN 
client  device,  in  range  of  the  access  point,  can  receive  radio  frequency  data  transmitted  to 
or  from  the  access  point.  Thus,  the  IEEE  802.11  specification,  which  applies  to  802.11a, 
802.11b  and  802.1  lg,  stipulates  some  security  requirements  as  described  in  the  next 
section. 

B.  WLAN  VULNERABILITY  AND  SECURITY 

Since  radio  waves  travel  through  ceilings,  floors,  walls,  windows  and  doors, 
unauthorized  persons  may  receive  the  transmitted  data  by  being  on  different  floors  or 
even  outside  the  physical  boundary  of  the  building.  If  a  WLAN  does  not  have  strict 
security  measures,  it  is  essentially  giving  wired-network  access  to  everyone  within  range, 
including  those  outside  of  the  building  [2],  Since  some  companies  deploy  WLANs 
without  all  security  features,  WLANs  attract  the  attention  of  hackers.  A  hacker,  with  a 
client  device  within  radio  frequency  range  of  an  access  point,  can  intercept  and  disrupt 
the  transmitted  data  if  the  data  is  not  encrypted. 

No  one  can  guarantee  a  completely  secure  WLAN  that  will  prevent  unauthorized 
users  from  penetrating  the  network.  This  does  not  mean  that  security  measures  are 
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pointless.  Ongoing  security  protections  must  be  established  and  frequently  updated.  If  the 
security  protection  is  deployed  and  left  static,  an  illegitimate  user  will  likely  get  access  to 
the  network  eventually. 

Two  main  security  approaches  for  WLANs  are: 

1 .  limit  the  signal  coverage  by  using  low  power  transmitters  for  a  predefined 
physical  zone  cleared  of  illegitimate  client  devices,  and 

2.  employ  authentication  and  privacy  features. 

The  first  method  can  be  achieved  by  either  measurement  or  simulation  using  wireless 
toolsets.  The  properties  of  the  antenna  and  the  environment  are  defined  in  the  software 
input  file  along  with  the  placement  of  antennas.  The  power  levels  for  coverage  of  an 
access  point  will  be  determined  after  running  enough  simulations.  By  antenna  selection 
and  placement,  the  signal  coverage  can  be  limited  to  the  area  occupied  by  legitimate 
users. 

The  second  approach  is  to  prevent  unintended  recipients  from  obtaining  access  to 
the  WLAN  when  they  are  within  the  coverage  of  the  network.  Authentication,  also  called 
access  control,  prevents  unauthorized  clients  from  sending  or  receiving  WLAN  data. 
WLAN  security  for  authentication  includes  the  combination  of 

•  Service  Set  Identifiers  (SSIDs), 

•  Wired  Equivalent  Privacy  (WEP)  keys, 

•  Media  Access  Control  (MAC)  authentication,  and 

•  open  or  shared  key  authentication. 

WLAN  privacy,  on  the  other  hand,  ensures  that  only  legitimate  users  can  use  the 
transmitted  data.  When  the  data  is  encrypted  with  a  key  that  can  be  used  only  by 
legitimate  recipients,  the  privacy  of  the  WLAN  is  considered  protected.  Encrypting 
WLAN  data  shows  that  the  integrity  of  the  data  remains  uncorrupted  during  the 
transmission  process  [3]. 

In  a  military  conflict,  it  might  be  necessary  to  “crash”  or  render  inoperable  an 
enemy’s  WLAN.  If  this  is  not  possible  by  hacking  into  the  network,  then  jamming  is 
another  alternative.  The  jammer  transmits  a  strong  noise  or  interference  signal  that 
prevents  the  WLAN  receiver  from  detecting  legitimate  network  signals. 
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To  be  effective,  the  jammer  must  transmit  relatively  high  power  or  be  very  close 
to  the  WLAN  access  point.  If  the  access  point  is  in  a  physically  secure  area,  then  high 
jammer  power  or  high  antenna  gain  is  required.  If  the  jammer  is  located  outside  of  the 
building,  then  the  jammer  propagation  losses  can  be  very  high.  On  the  other  hand,  if  there 
are  unprotected  window  openings,  then  the  jammer  signal  can  possibly  penetrate  into  the 
building  and  effectively  disrupt  the  network. 

C.  PREVIOUS  STUDY 

Indoor  propagation  was  previously  studied  by  Sumagaysay  [4]  and  Boukraa  [5]. 
The  Urbana  toolset  was  used  predict  the  contours  of  WLAN  power  levels  for  a  given 
building  model.  The  WLAN  signal  distribution  in  single  story  and  high-rise  buildings  at  a 
frequency  of  2.4  GHz  was  examined.  Boukraa  concluded  that  the  signal  levels  for  a  1- 
Watt  transmitter  could  only  be  detected  at  the  -70  dBm  level  within  two  floors  of  a  high- 
rise  structure.  Two  WLAN  vulnerabilities  were  apparent.  The  first  was  the  spread  of  the 
signal  outside  the  building  through  the  windows,  and  the  second  was  the  penetration  of 
the  signal  through  ceilings  and  floors. 

D.  OBJECTIVE  AND  APPROACH 

A  simple  jamming  transmitter  can  easily  interfere  with  WLANs  and  make 
communications  impossible.  Other  wireless  services,  such  as  Bluetooth,  that  operate  in 
the  same  frequency  range  can  also  interfere  with  WLAN  transmissions  and  reduce  the 
usable  bandwidth  of  the  WLAN  [3], 

The  purpose  of  this  thesis  is  to  examine  the  effects  of  noise  jamming  on  a  WLAN 
in  a  high-rise  building.  Specifically,  the  frequency  of  2.4  GHz  (IEEE  802.11b  and  IEEE 
802.1  lg  standard)  was  studied. 

It  was  assumed  that  the  security  policies  in  the  selected  high-rise  building  were 
very  stringent.  The  signal  coverage  was  limited  to  the  specific  zone  by  using  a  low  power 
transmitter  of  100  mW  [6],  which  was  located  in  the  most  interior  part  of  the  building.  It 
was  also  assumed  that  the  policy  did  not  allow  electronic  devices,  which  might  possibly 
interfere  with  the  WLAN.  If  the  entire  building  is  physically  secure,  then  the  closest  point 
for  jamming  the  WLAN  would  be  about  30  meters  away  from  the  building,  in  a  parking 
area.  Chapter  IV  describes  the  properties  of  the  jammer  in  detail. 
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E.  THESIS  OUTLINE 

Following  the  Introduction,  Chapter  II  discusses  radiowave  propagation  and 
interference  as  it  relates  to  wireless  propagation.  Chapter  III  gives  an  overview  of  the 
Rhino  software  which  is  used  to  create  the  building  model,  and  the  Urbana  software, 
which  is  used  to  perform  the  propagation  simulations.  Chapter  IV  describes  the  properties 
of  the  jammer.  Simulation  and  analysis  of  wireless  propagation,  interference,  and 
jamming  in  a  high-rise  building  is  also  contained  in  this  chapter.  Chapter  V  summarizes 
the  results  of  the  simulations  and  gives  recommendations  for  future  work. 
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II.  RADIO  WAVE  PROPAGATION 


Electromagnetic  waves  propagate  by  a  variety  of  mechanisms  that  are  generally 
categorized  as  reflection,  diffraction,  and  scattering.  In  an  indoor  environment,  where 
there  is  no  direct  line-of-sight  (LOS)  between  the  transmitter  and  the  receiver,  diffraction 
and  multiple  reflections  from  different  objects,  cause  the  electromagnetic  waves  to 
propagate  along  different  paths  and  have  different  arrival  times  and  phase  shifts. 
Multipath  fading  happens  because  of  the  interaction  between  these  waves.  [7] 

A.  PROPAGATION  MECHANISMS 

1.  Electromagnetic  Fields  and  Waves 

Electromagnetic  waves  consist  of  sinusoidal  electric  and  magnetic  fields,  which 
are  completely  described  by  Maxwell’s  equations: 


Y7  ?  dH 

VxE  = -u - 

dt 

(2.1) 

VxH  =  J  +  s  — 
dt 

(2.2) 

V-H  =  0 

(2.3) 

V-E  =  pv/s 

(2.4) 

where  E  (V/m)  is  the  electric  field  intensity,  H  (A/m) 

is  the  magnetic  field  intensity, 

J  (A/m2)  is  the  volume  current  density,  and  pv  (C/m3) 

is  the  volume  charge  density. 

Electrical  properties  of  a  medium  through  which  the  waves  are  propagating  are  specified 

by  its  constitutive  parameters,  such  as  permittivity  (s  in  F/m) ,  permeability  (//  in  H/m) , 

and  conductivity  (cr  in  S/m) . 

Assuming  a  e]COt  time  dependence,  the  complex  permittivity  and  permeability  are 

e  =  e'-je" 

(2.5) 

p  =  p'-jp" . 

(2.6) 

If  the  medium  is  lossy,  then  the  imaginary  terms  are  nonzero. 
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The  relative  values  for  s  and  //  are 


*  =  e0er 

(2.7) 

M  =  Mo  Mr 

(2.8) 

where  er  is  the  relative  permittivity,  and  /./.  is  the  relative  permeability.  Both  can  be 

complex  in  accordance  with  Equations  (2.5)  and  (2.6).  In  free  space,  which  is  a  lossless 
medium 

s  =  e0  (2.9) 

where  e0  =  8.8542 xlCT12  F/m,  and 


M  =  Mo  (2.10) 

where  ju0  =4^x10  7  H/m.  The  majority  of  the  materials  in  an  indoor  environment  are 
nonmagnetic.  Thus,  ju  =  /u0  is  assumed  for  further  calculations. 

In  terms  of  the  conductivity  of  the  medium,  the  complex  relative  dielectric 
constant  of  the  material  is  [8] 


£•„  =  £„ 


cosn 


(2.11) 


where  a  is  the  conductivity  of  the  material  (S/m),  co  =  2nf  is  the  radian  frequency 
(rad/sec),  and  /  is  the  frequency  (Hz).  Finally,  the  intrinsic  impedance  of  the  medium  is 


71  = 


(2.12) 


The  wave  equations  can  be  derived  for  the  electric  and  magnetic  fields.  They 
solve  Maxwell’s  equations  for  a  source  free  region 


V2E-y2E  =  0 

(2.13) 

V2H-y2H  =  0. 

(2.14) 
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Plane  waves  are  the  simplest  solutions  to  the  wave  equations.  A  plane  wave  propagating 
in  the  z  direction  is  given  by  the  phasor  expression 


E(z)  =  xE0e~rt 


(2.15) 


where  E0  is  the  electric  field  at  the  source,  and  y  is  the  propagation  constant,  which  is 

y  =  j(oJJie=a  +  jj3  (2.16) 

2n 

where  a  is  the  attenuation  constant  (Np/m),  p  =  —  is  the  phase  constant  (rad/m),  and 

/l  is  the  wavelength.  The  phase  velocity  of  electromagnetic  wave  in  a  medium  depends 
on  the  permittivity  and  permeability  of  the  medium 


co 


(2.17) 


In  free  space,  up=c=  ,  =  2.998 x  108  m/s . 

"V  //Ao 

2.  Reflection  and  Refraction 

When  an  electromagnetic  wave  propagates  in  one  medium  and  impinges  upon 
another  medium  with  different  electrical  properties,  reflection  occurs.  When  the  boundary 
has  large  dimensions  relative  to  the  wavelength  of  the  propagating  wave,  then  high 
frequency  approximation  (i.e.,  ray  tracing)  applies. 

The  wave  will  also  be  partially  refracted  into  the  second  medium  unless  it  is  a 
perfect  conductor.  If  the  second  medium  is  a  perfect  electric  or  magnetic  conductor,  there 
will  be  no  refraction,  and  the  entire  incident  wave  is  reflected  back  into  the  first  medium. 
The  reflected  and  refracted  wave  intensities  are  also  related  to  the  properties  of  the 
incident  wave,  such  as  the  wave  polarization,  angle  of  incidence,  and  the  frequency  of  the 
propagating  wave.  [7]  A  sample  ray  diagram  of  a  plane  wave  incident  on  an  interface 
between  two  dielectrics  is  illustrated  in  Figure  1 . 
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Normal 


Figure  1 .  A  sample  of  reflected  and  refracted  waves 

The  angle  of  the  reflected  wave  is  equal  to  the  angle  of  incidence  by  Snell’s  law 
of  reflection 

0t=0r.  (2.18) 

This  is  referred  to  as  specular  reflection.  The  index  of  refraction  of  a  medium  is  the  ratio 
of  the  phase  velocity  in  free  space  to  the  phase  velocity  in  the  medium.  By  Snell’s  law  of 
refraction 

nx  sin  6t  =  n2  sin  6t  (2.1 9) 

where  nx  and  n2  are  the  indices  of  refraction  of  media  1  and  2,  respectively.  The  index  of 

refraction  is  related  to  the  dielectric  constant,  n  =  ^£~r  .  When  «,  <n2,  it  follows  that 

sin  <9(  >  sin  6t  and  also  that  dt  >  6t .  Thus,  it  can  be  said  that  the  ray  entering  a  higher 

index  medium  bends  toward  the  normal  as  in  Figure  2. (a).  The  reverse  is  also  true,  that  is, 
the  ray  entering  a  lower  index  medium  bends  away  from  the  normal  as  in  Figure  2.(b)  [9], 
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(a) 


(b) 


Figure  2.  Bending  of  rays  at  an  interface,  (a)  the  ray  entering  a  higher-index 
medium,  (b)  the  ray  entering  a  lower  index-medium 

3.  Diffraction 

Diffraction  occurs  when  there  is  a  surface  that  has  a  discontinuity  between  the 
transmitter  (Tx)  and  the  receiver  (Rx),  and  obstructs  the  direct  path  of  the 
electromagnetic  wave.  The  secondary  wavelets  created  at  the  edge  propagate  behind  the 
obstacle  into  the  shadow.  Thus,  there  can  be  sufficient  field  strength  for  reception  in  the 
shadowed  region. 

The  intensity  of  diffraction  depends  on  the  geometry  of  the  obstacle,  and 
amplitude  and  polarization  of  the  incident  electromagnetic  wave.  A  sample  of  edge 
diffraction  is  illustrated  in  Figure  3. 
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Figure  3.  Geometry  of  wedge  diffraction 

4.  Scattering 

When  electromagnetic  waves  arrive  at  the  receiver,  the  received  signal  is  often 
stronger  than  what  is  expected  based  on  specular  reflection  and  diffraction.  Scattering 
occurs  when  a  electromagnetic  wave  impinges  upon  a  rough  surface  or  object  which  has 
small  dimensions  relative  to  the  wavelength.  This  is  because  scattered  waves  propagate  in 
many  directions  due  to  isotropic  scattering.  In  an  indoor  environment,  such  as  in  offices, 
a  variety  of  small  objects  (e.g.,  cups,  computer  mouse,  staplers,  etc.)  scatter  waves  in  all 
directions.  [7] 

B.  ANTENNA  FUNDAMENTALS 

An  antenna  may  be  defined  as  a  transducer  that  provides  a  means  for  radiating  or 
receiving  electromagnetic  waves.  Antennas  are  made  in  different  sizes  and  shapes,  and 
are  used  in  various  areas,  such  as  [10] 

•  radiowave  communication  systems  (e.g.,  WLANs), 

•  cellular  systems, 

•  radio  and  television  broadcasting  and  reception,  and 

•  radar  systems. 

An  isotropic  antenna  is  a  hypothetical  antenna  that  has  a  uniform  radiation  pattern  in  all 
directions.  It  is  used  as  a  reference  for  defining  the  directive  properties  of  all  antennas. 
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1.  Antenna  Radiation  Pattern 

The  antenna  radiation  pattern  is  the  angular  variation  of  radiation  fields  at  a  fixed 
distance  from  the  transmitting  antenna  [11],  The  antenna  pattern  is  represented  by 
F  (<9, (f)  for  the  field  and  P{0,(/))  for  the  power.  The  general  normalized  power  pattern  in 
decibels  is 

P(0,tf>)  =  201og  [dB]  (2.20) 

■^max 

where  Fmax  is  the  main  lobe  magnitude. 

The  antenna  pattern,  as  shown  in  Figure  4,  has  main,  side  and  back  lobes.  Side 
and  back  lobes  together  are  also  called  minor  lobes.  The  side  lobe  level  (SLL)  is  the  ratio 
of  the  side  lobe  peak  to  the  main  lobe  peak 

SLL  =  201og^  [dB]  (2.21) 

F 

max 

where  FSLL  is  the  highest  side  lobe  magnitude.  The  angular  width  of  the  main  lobe  is  the 

angle  between  first  nulls  and  is  generally  represented  by  half-power  beam  width  (HP  or 
HPBW).  The  HPBW  is  also  known  as  the  3-dB  beam  width, 

HpBW=  (2.22) 

A  sample  of  an  antenna  radiation  pattern  is  illustrated  in  Figure  4. 
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Figure  4.  Sample  polar  plot  of  antenna  radiation  pattern  (From  Ref.  [1 1].) 


2.  Antenna  Directivity 

The  directivity  is  the  ratio  of  maximum  normalized  power  density  in  the  main 
beam  direction  of  the  radiation  pattern  to  the  average  power  density  of  the  antenna 


(2.23) 


where  S  (W/m2)  is  the  maximum  power  density,  and  Sj  is  the  average  power  density.  It 

is  assumed  that  all  of  the  power  in  an  antenna  is  radiated.  When  describing  the  pattern 
properties  of  an  antenna,  the  radiation  pattern  is  compared  to  spatial  distribution  of  an 
isotropic  antenna.  For  an  isotropic  antenna,  the  directivity  is  equal  to  one.  A  sample  of  an 
antenna  radiation  pattern  and  its  directivity  is  illustrated  in  Figure  5. 


12 


Pattern  F(9 ,  #)  of 
a  real  antenna 


Isotropic 

pattern 


Figure  5.  Illustration  of  radiation  pattern  and  directivity  (From  Ref.  [11].) 

3.  Antenna  Gain 

The  directivity  is  reduced  due  to  losses  in  the  antenna.  This  loss  defines  the 
efficiency  of  the  antenna,  e.  The  gain  is  then  the  product  of  the  directivity  and  the 
efficiency.  If  the  efficiency  of  the  antenna  is  100%,  the  gain  is  equal  to  directivity.  In 
practice,  the  gain  of  a  lossless  antenna  is 

A  jr  A 

G  =  (for4»T2)  (2.24) 


where  Ae  is  the  effective  area  of  the  antenna.  The  gain  of  an  antenna  is  also 
approximately  equal  to  [12] 


G 


26,000 

- e 


(2.25) 


where  0B  and  4  are  the  elevation  and  azimuth  HPBWs,  respectively,  in  degrees. 

4.  Antenna  Polarization 

The  antenna  polarization  refers  to  the  vector  direction  of  the  wave  radiated  by  the 
transmitting  antenna  in  terms  of  the  time-varying  orientation  of  the  electric  field.  The  tip 
of  the  electric  field  vector  traces  a  path  in  space,  and  generally,  it  is  elliptical  [10].  This  is 
called  elliptical  polarization.  If  the  path  is  a  circle  or  a  straight  line,  then  the  polarization 
is  called  circular  or  linear  polarization,  respectively. 
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In  most  wireless  applications,  the  transmitting  and  the  receiving  antennas  are 
linearly  polarized.  Polarization  mismatch  loss  occurs  if  the  transmitting  and  receiving 
antennas  are  not  similarly  polarized. 

C.  FREE-SPACE  LINK  EQUATION 

Assuming  no  polarization  mismatch  loss  and  no  obstruction  between  the 
transmitting  and  the  receiving  antennas,  as  shown  in  Figure  6,  the  received  signal  power 
in  free-space  can  be  expressed  by  Friis’s  transmission  formula  [11] 


P  = 


(4 


(2.26) 


where  Pt  is  the  transmitted  power  in  watts,  d  is  the  distance  between  the  transmitting 
and  the  receiving  antennas,  and  Gt  and  Gr  are  the  gains  of  transmitting  and  the  receiving 
antennas,  respectively. 


Tx 


Pr,Gr 


Figure  6.  A  free-space  communication  link 

The  effective  isotropically  radiated  power  (EIRP)  is  the  transmitted  power  into 
the  antenna  multiplied  by  the  power  gain  of  the  transmitting  antenna  in  a  given  direction 
relative  to  an  isotropic  antenna 

EIRP  =  PtGt  .  (2.27) 

It  is  also  called  the  effective  radiated  power  (ERP).  Usually,  the  free-space  link  equation 
is  written  in  decibels 

Pr=P,+G,  +  Gr-(L„  +  L)  [dB]  (2.28) 
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where  Lfs  and  L  are  the  free-space  path  loss  and  the  other  losses,  respectively.  The  free- 

space  path  loss  is  proportional  to  the  inverse  square  of  the  distance  between  transmitting 
and  receiving  antennas  and  is  frequency  dependent 


Lfs  =10  log 


f  And2  ^ 
V  ?  7 


[dB], 


(2.29) 


D.  INDOOR  PROPAGATION 

Indoor  propagation  is  affected  by  reflection,  diffraction,  and  scattering  as  it  is  for 
outdoor  propagation.  However,  characterizing  indoor  propagation  differs  from  outdoor 
propagation  modeling  since  the  distances  are  much  smaller  and  the  environmental 
variability  is  much  greater.  The  building  type,  construction  materials,  and  indoor  layout 
(open  or  closed  doors)  contribute  to  the  received  signal  level  and  indoor  path  loss  [7]. 

1.  Log-Distance  Path  Loss  Model 

For  both  outdoor  and  indoor  propagation,  the  average  received  signal  power 
decreases  logarithmically  with  distance  [7] 


PL  =  PL(d0  )  +  lOn  log 


f d_^ 

V^o  J 


[dB] 


(2.30) 


where  PL  is  the  path  loss,  d0  is  the  reference  distance,  d  is  the  distance  between  the 

transmitter  and  the  receiver,  and  n  is  the  path  loss  exponent  which  depends  on  the 
specific  propagation  environment.  The  exponent  n  will  have  a  larger  value  when 
obstructions  are  present  between  the  transmitter  and  the  receiver.  The  bars  in  Equation 
(2.30)  denote  the  average  path  loss  values  for  a  given  value  of  d  .  Table  1  shows  the  path 
loss  exponents  (n)  for  different  environments. 


Environment 

Path  Loss  Exponent 

Free  space 

2 

Urban  area  cellular  radio 

2.7  to  3.5 

Shadowed  urban  cellular  radio 

3  to  5 

In  building  line-of-sight  (LOS) 

1.6  to  1.8 

Obstructed  in  building 

4  to  6 

Obstructed  in  factories 

2  to  3 

Table  1 .  Path  loss  exponents,  n ,  for  different  environments  (From  Ref.  [7].) 
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2.  Attenuation  Model 

The  building  type  and  obstacles  in  the  building  degrade  the  performance  of 
wireless  systems  indoors  [7],  A  common  attenuation  factor  model  is 


PL(d)  =  PL(d0)  +  20 log  —  +ad  +  FAF  +  YJPAF  [dB] 


V^o  J 


(2.31) 


where  a  is  the  attenuation  constant  in  dB/m ,  FAF  is  the  floor  attenuation  factor,  PAF 
is  the  partition  attenuation  factor  for  an  obstruction  which  blocks  the  LOS  between  the 
transmitter  and  the  receiver.  As  shown  in  Equation  (2.31),  in  multi-floor  buildings,  the 
path  loss  depends  on  the  ffee-space  loss  plus  an  additional  loss,  which  increases 
exponentially  with  distance. 

3.  Ray  Tracing  Model 

New  methods  for  predicting  signal  coverage  in  indoor  environments  involve  some 
computational  electromagnetic  tools,  such  as  Urbana,  which  is  described  in  detail  in  the 
next  chapter.  This  software  uses  a  three-dimensional  (3-D)  ray- tracing  engine,  which  also 
includes 

•  antenna  pattern  and  polarization  effects, 

•  fading  and  co-channel  interference,  and 

•  multipath  and  diffraction  effects. 

With  the  help  of  these  computer-aided  design  tools,  predictions  and  measurements  of  the 
signal  coverage  may  be  simultaneously  viewed,  manipulated  and  optimized  [7]. 

E.  INTERFERENCE  AND  COMMUNICATION  JAMMING 

Interference  occurs  when  two  signals  at  the  same  frequency  arrive  at  a  single 
receiver.  As  shown  in  Figure  7,  one  of  these  signals  is  considered  the  desired  signal  and 
the  other  is  an  interfering  signal.  The  interfering  signal  could  be  either  intentional  or 
unintentional.  If  it  is  intentional,  then  the  objective  of  the  interfering  signal  is  to  reduce 
the  quality  of  the  desired  signal  to  an  unusable  level  [13]. 
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Figure  7.  Sample  link  from  each  transmitter  to  the  receiver  (From  Ref.  [13].) 


The  power  difference  between  the  desired  signal  and  the  interfering  signal, 
assuming  that  the  receiving  antenna  has  the  same  gain  to  both  signals,  is  written  in 
decibels 


where 


Ps-Pj=  ERPS  -  ERPI  +  20  log  (d, )  -  20  log  ( ds  )  [dB] 


(2.32) 


Ps  is  the  received  power  from  the  desired  source, 

Pj  is  the  power  from  the  interfering  source, 

ERPS  is  the  effective  radiated  power  of  the  desired  signal, 

ERPf  is  the  effective  radiated  power  of  the  interfering  signal, 
ds  is  the  path  distance  to  the  desired  antenna,  and 
dj  is  the  distance  to  the  interfering  antenna. 

The  same  link  equation  can  be  applied  to  the  jammer  as  for  the  desired 

transmitter.  The  ratio  of  received  power  gives  the  signal-to-jam  ratio  i^/jj ,  which  is 
then  [14] 


S//J  =  ERPs-ERPJ+GR-GRJ+20log(dJ)-20\og(ds )  [dB]  (2.33) 


where  GR  and  GRJ  are  the  gains  of  the  receiving  antenna  toward  the  desired  transmitter 
and  toward  the  jammer,  respectively.  If  the  gain  of  the  receiving  antenna  is  omni¬ 
directional  in  azimuth,  GR  and  GRJ  cancel  each  other.  Otherwise,  if  the  receiving  antenna 
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has  a  narrow  beam  and  is  pointed  at  the  desired  transmitter,  then  the  jammer  has  to 
transmit  enough  power  to  overcome  the  difference  in  the  antenna  gain. 

In  this  thesis,  the  objective  is  to  inject  an  interference  signal,  as  in  noise  jamming, 
into  the  receiver  such  that  the  actual  signal  is  submerged  in  the  interference  [15].  This 
type  of  jamming  is  also  known  as  denial  jamming. 

This  chapter  is  mainly  concerned  radiowave  propagation  mechanisms,  antenna 
fundamentals,  the  ffee-space  link  equation,  and  interference  and  communication 
jamming.  The  next  chapter  discusses  the  software  tools,  which  are  used  for  the  model  and 
the  simulations. 
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III.  SOFTWARE  DESCRIPTIONS 


A.  INTRODUCTION 

Urbana  is  a  Unix-based  program  developed  by  Science  Application  International 
Corporation  (SAIC).  The  main  components  of  Urbana  are  Cifer,  XCell,  and  Urbana 
[16]. 

Urbana,  being  a  powerful  Computational  Electromagnetic  (CEM)  tool, 
determines  the  signal  level  for  wireless  networks  with  predefined  inputs,  such  as, 
building  geometry,  antenna  types,  frequency,  and  polarization.  There  is  capability  to 
analyze  [17] 

•  urban  propagation, 

•  indoor  and  outdoor  propagation, 

•  collision  avoidance  radar,  and 

•  antenna  studies. 

The  purpose  of  this  research  is  to  analyze  the  effects  of  noise  jamming  on  wireless 
systems  operating  inside  a  high-rise  building. 

B.  SOFTWARE  TOOLS 

1.  Rhino 

Rhino  Version  2.0  software,  which  runs  on  Windows  platforms,  was  chosen  to 
generate  the  high-rise  building.  It  is  a  very  powerful,  and  user  friendly  CAD  software. 
Rhino  is  a  very  effective  toolset  to  create,  analyze,  render,  edit,  and  build  a  3-D  model 
precisely.  [18]  The  individual  components  of  the  model  for  this  project  were  created 
separately  by  Rhino,  such  as,  the  outer  walls,  the  inner  walls,  the  doors,  the  windows, 
and  the  elevators.  This  allowed  the  materials  for  each  group  to  be  changed  independent  of 
the  others. 

Rhino  supports  many  file  formats.  For  this  project,  the  * .iges  format  was  chosen. 
This  format  was  readable  by  Cifer,  as  described  in  the  next  section.  A  sample  of  the  four- 
view  Rhino  editor  is  illustrated  in  Figure  8. 
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Figure  8.  A  sample  Rhino  GUI  editor 

2.  Cifer 

Cifer  is  primarily  used  to  convert  files  imported  from  CAD  software  programs 
into  the  facet  file  format  required  by  Urbana.  All  the  predefined  *.iges  files  for  the 
individual  building  components  imported  from  Rhino  were  converted  to  *. facet  files  with 
different  ICO  AT  numbers.  Thus,  each  ICO  AT  number  corresponds  to  a  unique  material. 
Zero  is  the  ICOAT  value  for  a  perfect  electric  conductor  (PEC). 

The  model  is  completed  by  combining  all  of  the  individual  facet  files  into  one. 
For  instance,  the  combination  of  the  outer  walls  and  the  inner  walls  is  subsequently 
combined  with  the  doors  and  so  forth,  to  obtain  a  single  file  for  the  entire  building. 

Cifer  is  also  capable  of  extracting  edges  from  a  facet  file.  The  consolidate  mode 
is  also  used  to  consolidate  facets  and  edges.  A  sample  graphical  user  interface  (GUI)  is 
illustrated  in  Figure  9. 
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Figure  9.  A  sample  Cifer  GUI  editor 

3.  XCell 

XCell  is  used  to  view  and  modify  the  facet  files  created  by  Cifer  and  to  view  the 
signal  contours  generated  by  the  programs  f2f  or  f2fd.  One  can  view  the  3-D  perspective 
of  signal  strengths  in  and  around  the  building  model  in  color-coded  plots.  It  is  also 
possible  to  view  facet  files  together  with  antenna  files  and  edge  files  at  the  same  time.  A 
sample  GUI  for  XCell  is  illustrated  in  Figure  10. 
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Figure  10.  A  sample  XCell  GUI  editor 

4.  Urbana 

Urbana  uses  a  3-D  ray  tracing  engine  to  characterize  electromagnetic  wave 
propagation.  The  ray  tracing  engine  has  a  high-speed  algorithm  to  implement  Geometric 
Optics  (GO),  Physical  Optics  (PO)  and  the  Geometrical  Theory  of  Diffraction  (GTD).  A 
sample  Urbana  file  appears  in  the  Appendix.  The  key  inputs  to  the  Urbana  code 
(* .urjnput)  are 

•  facet  model  for  the  CAD  geometry  file  (*. facet), 

•  frequency, 

•  antenna  file,  which  includes  antenna  pattern,  power,  and  polarization 
(*  .antenna), 

•  observation  point  file  (*.txt), 

•  edge  file  for  edge  diffraction  (*  .edge),  and 

•  building  material  properties,  such  as,  concrete,  ground,  wood,  glass,  and 

dielectrics. 

The  key  output  to  the  Urbana  code  is  a  composite  field  file  (*. field)  for  the 

observation  plane,  which  is  defined  in  the  input  file. 

22 


a.  Generating  Input  Data  Files  (Step  1) 

The  building  model  is  represented  by  surfaces  of  triangular  facets.  As 
described  in  Section  2,  Rhino  is  used  to  generate  the  high-rise  building  model  in  *.iges 
format.  Cifer  is  then  used  to  convert  the  output  file  to  the  DEMACO  facet  file  format. 
Each  facet  file  is  characterized  by  the  coating  material  defined  by  the  input  variable 
ICOAT.  Cifer  is  also  used  to  extract  edges  (*.edge)  from  the  facet  file.  Urbana  also 
requires  a  set  of  observation  points  at  which  the  field  is  to  be  computed.  To  build  an 
observation  point  list,  a  Matlab  script  file  is  run  and  the  output  file  is  saved  in  ASCII 
*.txt  format.  Step  1  is  illustrated  in  Figure  11. 


Generate  building  structures 

Modify  or  combine  structures; 
add  or  change  material;  translate, 
rotate  or  scale  dimensions 

Interactively  and  graphically 
modify  structures  and  materials 

Build  observation  point  list 
and  facet  file  for  viewing 


(obv).list,  (obv). facet 


Figure  11.  Generating  input  data  file  for  Urbana  (From  Ref.  [19].) 

b.  Generate  Urbana  Input  File  and  Run  Urbana  (Step  2) 

The  Urbana  input  file  (* .urjnpui)  is  composed  of  the  key  input 
parameters  and  the  file  names  described  above.  After  executing  Urbana,  the  output  is 
summarized  in  the  composite  field  file  (*. field).  This  step  is  illustrated  in  Figure  12. 
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Urbana  input  file  includes  the 
file  names  generated  in  STEP  1 ; 
defines  link  parameters,  material 
electrical  parameters,  and  EM 
propagation  parameters 


(geo),  facet 
(runname). urinput  (geo). edge 

(obv).list 


urbana.x 


I 

(runname).  field 


Figure  12.  Process  for  running  Urbana  (From  Ref.  [19].) 
c.  Post  Processing  (Step  3) 

After  getting  the  *. field,  file,  the  command  f2f  is  used  to  convert  it  to  a 
*. facet  file  for  viewing  the  signal  back  in  XCell.  There  are  some  selections  with  this 
command,  such  as,  magnitude  of  electric  field  level,  antenna  power  level,  color  bar 
settings,  clip  minimum-maximum,  and  range  minimum-maximum.  XCell  is  used  to  view 
the  *. facet  file.  This  step  is  illustrated  in  Figure  13. 


Convert  Urbana  output  field  values 
to  viewable  (facet)  contours 


View  signal  contours  and 
geometry  structures  simultaneously 


(runname).  field 

i 

f2f.x 

|  (runname). facet 

(obv).  facet 
(geo),  facet 
(xxx).xedge_pref 


xcell.x 


Figure  13.  Urbana  post-processing  (From  Ref.  [19].) 

C.  DESIGNING  MODEL 

As  discussed  in  Section  B,  the  model  is  built  using  the  Rhino  software.  The 
building  is  a  ten-story  office  or  hotel  building.  The  first  floor  has  a  larger  area  than  the 
other  floors.  The  dimensions  of  first  floor  are  illustrated  in  Figure  14. 
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x  y 


Figure  14.  First  floor  lobby  area 


Floors  2  through  9  are  identical.  They  consist  of  two  rows  of  rooms  (or  offices) 
with  elevators  and  stairwell.  Also,  restrooms  are  located  at  one  end  of  the  building,  as 
shown  in  Figure  15. 


Figure  15.  Floors  two  through  nine  (office  areas  or  hotel  rooms) 

The  roof  has  an  open  patio  area  and  an  enclosed  restaurant  area.  The  dimensions 
of  the  roof  are  illustrated  in  Figure  16. 
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Figure  16.  Tenth  floor  with  patio  and  restaurant 

Finally,  a  flat  ground  plane  is  located  under  the  building  to  simulate  the  ground 
with  the  dimensions  as  shown  in  Figure  17. 


Figure  17.  The  whole  building  model  with  the  ground 

The  model  is  composed  of  six  different  materials,  and  each  one  is  given  a 
different  ICOAT  value  after  it  is  imported  into  Cifer.  All  parts  are  then  combined  by 
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using  the  “combine  two  facet  files”  command  in  Cifer.  The  building  consists  of  the 
following  materials: 

•  The  outer  walls  (ICOAT=l)  which  are  concrete  and  0.3  m  thick. 

•  The  inner  walls  (ICOAT=2)  also  concrete  but  0.1  m  thick. 

•  The  doors  (ICOAT=3)  which  are  wood  and  0.05  m  thick. 

•  The  windows  (ICOAT=4)  are  glass  and  0.025  m  thick. 

•  The  elevator  doors  (ICOAT=0)  are  PEC  and  0.05  m  thick. 

•  The  ground  (ICOAT=6)  is  considered  as  a  0.1  m  thick  concrete  layer  over 
a  semi-infinite  conductor. 

Each  component  has  its  own  permeability,  permittivity,  resistivity,  and  thickness 
values  specified  in  the  *.ur_input  file  as  shown  in  Table  2. 


COATING 

MATERIAL 

Material 

ICOAT 

Permeability 

(A',  A") 

Permittivity 

(e',en) 

Resistivity 

(ohm) 

Thickness 

(meters) 

Outer  Walls 

Concrete 

1 

1,0 

10.1, -0.5 

le30 

0.3 

Inner  Walls 

Concrete 

2 

1,0 

10.1, -0.5 

le30 

0.1 

Doors 

Wood 

3 

1,0 

3,  -0.67 

le30 

0.05 

Windows 

Glass 

4 

1,0 

6,0 

le30 

0.025 

Elevator  Doors 

Metal 

0(PEC) 

1,0 

1,0 

0 

0.05 

Ground 

Concrete 

6 

1,0 

3,0 

le30 

0.1 

Table  2.  Properties  of  each  coating  material 


This  chapter  dealt  mainly  with  explaining  different  software  tools  used  to  perform 
the  simulations.  The  next  chapter  concerns  creating  the  jammer  antenna  and  seeing  the 
effects  of  noise  jamming  on  the  model. 
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IV.  URBANA  SIMULATIONS,  ANALYSIS  AND  RESULTS 


The  previous  chapter  described  the  software  used  to  develop  the  high-rise 
building  model  and  a  brief  introduction  to  propagation  modeling.  This  chapter  describes 
the  WLAN  and  the  jammer  system  and  antenna  properties  as  well  as  the  placements,  the 
simulations  and  their  results.  Indoor  propagation  and  signal-to-jam  ratio  (SJR)  are  also 
evaluated. 

The  first  step  is  to  study  the  indoor  propagation  and  signal  distribution  for  the 
entire  building  model  without  the  jammer.  Next,  the  effects  of  the  jammer,  which  could 
possibly  have  different  HPBWs  and  power  levels,  are  calculated  inside  the  building 
model.  Finally,  the  presence  of  the  jammer  on  the  WLAN  system  is  studied.  By 
examining  the  SJR,  the  jamming  contours  are  generated  for  the  WLAN  network. 

A.  WLAN  TRANSMITTER  ANTENNA 

Currently,  most  WLAN  devices  operate  with  vertically  polarized  antennas.  A 
vertically  polarized  antenna  generally  has  a  slightly  higher  signal  level  than  a  horizontally 
polarized  antenna  in  an  indoor  environment  [20],  In  the  simulations,  the  transmitter 
antenna  is  considered  as  an  access  point  and  is  located  in  the  most  interior  part  of  the 
building  model,  i.e.,  in  the  middle  of  the  fifth  floor.  Also,  to  provide  network  coverage 
where  needed,  an  omni-directional  WLAN  antenna  is  used.  A  transmitter  with  a  power  of 
100  mW  with  802.1  lb/g  standard  has  a  data  rate  from  11  Mbps  to  54  Mbps.  The 
maximum  range  of  a  typical  access  point  configured  with  a  2.2  dBi  dipole  antenna  for 
802.1  lb/g  is  shown  in  Table  3  [21].  Thus,  2.4  GHz,  an  omni-directional  half-wave  dipole 
antenna,  transmitting  a  power  of  100  mW  is  chosen  for  the  simulations.  Figure  18  shows 
the  3-D  radiation  pattern  of  a  dipole  antenna. 


29 


Typical  Indoor  Data  Rate 

Range 

54  Mbps 

27  m  (90  ft) 

18  Mbps 

54  m  (180  ft) 

1 1  Mbps 

48  m  (150  ft) 

6  Mbps 

91  m  (300  ft) 

1  Mbps 

124  m  (410  ft) 

Table  3. 


The  maximum  range  and  the  data  rate  of  a  typical  access  point  (From  Ref. 

[21].) 


V=|E|Sir(t)sm(pJ  '1  -1  X=|E|Sin(l]CO$[p) 

Figure  18.  Normalized  radiation  pattern  of  a  dipole  antenna  (From  Ref.  [19].) 

B.  JAMMER  ANTENNA 

Since  the  entire  building  is  physically  secure  against  possible  inside  jammers  as 
mentioned  in  Chapter  I,  the  jammer  is  considered  to  be  located  about  30  meters  away 
from  the  building,  in  the  parking  area.  The  effects  of  three  different  jammers  each  with 
two  different  power  levels  are  studied  against  the  building  model.  Jammer  1,  jammer  2 
and  jammer  3  have  10°-HPBW,  20°-HPBW,  and  30°-HPBW,  respectively.  Each  jammer 
antenna  transmits  10  W  and  100  W. 

It  is  also  considered  to  use  a  vertically  polarized  directional  antenna  for  the 
jammer  to  obtain  maximum  jamming  effectiveness  against  the  WLAN  in  the  building.  A 
Matlab  code  that  plots  the  pattern  of  the  jamming  antenna  was  developed.  The  antenna 
radiation  pattern  of  the  jammer  is  zero  in  the  rear  hemisphere  of  the  jammer  antenna.  The 
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user  simply  inputs  the  HPBW,  and  the  generated  ASCII  pattern  file  is  exported  to 
Urbana  in  *  dat  file  format.  The  Matlab  code  appears  in  the  Appendix. 

Jammer  1  antenna  has  a  0.3650  meter  radius  and  25.2713  dBi  of  directivity. 
Sample  XCell  views  of  the  10° -HPBW  jammer  antenna  pattern,  with  10  W  and  100  W, 
are  illustrated  in  Figure  19. 


(90m,  -70m)  (-30m,  -70m) 


Figure  19.  The  radiation  pattern  of  10°  -HPBW  jamming  antenna  with  (a)  10  W,  (b) 

100  W 

Jammer  2  antenna  has  a  0.1825  meter  radius  and  19.2507  dBi  of  directivity. 
Sample  XCell  views  of  the  20°  -HPBW  jammer  antenna  pattern,  with  10  W  and  100  W, 
are  illustrated  in  Figure  20. 
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(90m,  -70m)  (-30m,  -70m) 


+  20  dBm 


■ 

I 


-  90  dBm 


Figure  20.  The  radiation  pattern  of  20°  -HPBW  jamming  antenna  with  (a)  10  W,  (b) 

100  W 


Jammer  3  antenna  has  a  0.1217  meter  radius  and  15.7288  dBi  of  directivity. 
Sample  XCell  views  of  the  30° -HPBW  jammer  antenna  pattern,  with  10  W  and  100  W, 
are  illustrated  in  Figure  21. 


Figure  2 1 .  The  radiation  pattern  of  30°  -HPBW  jamming  antenna  with  (a)  1 0  W,  (b) 

100  W 
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For  the  simulations,  each  jammer  antenna  is  then  tilted  30°  with  respect  to  the 
ground  to  aim  the  main  beam  approximately  at  the  middle  of  the  building.  Figure  22 
shows  the  10°  -HPBW  jammer  antenna  pattern  when  tilted  up  to  the  center  of  the  building 
model. 


+  20  dBm 


Figure  22.  The  plot  of  tilted  jammer  antenna  radiation  pattern 

C.  OBSERVATION  PLANE 

The  previous  chapter  explains  that  Urbana  requires  a  set  of  observation  points  at 
which  the  electric  field  is  to  be  computed.  The  observation  point  file  can  be  created  either 
by  bobv.x  or  a  Matlab  code.  For  this  research,  a  Matlab  script  file  was  run  and  the 
output  file  was  saved  in  ASCII  *.txt  format.  Then,  the  *.  txt  file  was  exported  to  Urbana. 
This  process  was  done  for  the  first  floor  through  the  tenth  floor,  and  observation  planes 
were  created  1  meter  above  each  floor  level.  For  each  observation  plane,  a  step  of  0.5 
meter  between  points  was  used.  The  dimensions  of  each  observation  plane  are  75  meters 
by  40  meters.  The  Matlab  code  to  generate  the  observation  plane  appears  in  the 
Appendix. 

D.  URBANA  SIMULATIONS 

It  obvious  that  the  indoor  environment  is  very  dynamic  when  WLANs  are 
considered.  Especially,  people  and  their  movements,  the  movement  of  objects,  the  usage 
of  personal  electronic  devices  that  interfere  with  the  WLAN,  and  the  opening  and  closing 
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windows  and  doors  will  affect  the  signal  distribution  in  an  indoor  environment.  For  this 
reason,  a  fixed  picture  of  the  building  model  is  used  to  develop  an  approximation  of  the 
signal  distribution. 

For  the  simulations,  a  high-fidelity  model  with  closed  windows  is  studiedThe 
internal  doors  are  considered  open  for  all  cases. 

1.  High-Fidelity  Model 

First,  the  building  model  is  considered  with  the  realistic  materials  described  in 
Table  2.  The  power  of  the  transmitting  antenna  is  a  fixed  value  of  100  mW.  The  input 
parameters  and  the  values  for  the  Urbana  input  file  are  shown  in  Table  4. 


Input  Parameters 

Response  Values 

Input  file 

model#real.ur_input 

Facet  file 

modelfinal.facet 

Edge  file 

model0136.edge 

Observation  plane  file 

obv#.txt 

Antenna  type 

XU  omni-directional  dipole 

Antenna  power  (W) 

1  (linear  value  for  calculation) 

Antenna  polarization 

vertical 

Antenna  frequency 

2.4  GHz 

Antenna  coordinates  (x,y,z) 

30.0,  15.0,  16.9  (meters) 

Propagation  mechanism 

GO 

Edge  diffraction 

Yes  (UTD) 

Ray  angular  intervals  (degree) 

2 

Max  ray  bounces 

10 

Coating  materials 

ICOAT=  0,  1,  2,  4,  6  (as  described  in  Table  2) 

Table  4.  Input  parameters  for  the  high-fidelity  model  with  closed  windows 


After  running  Urbana  for  the  each  observation  plane,  the  model#real. field  file  is 
generated  for  the  post  processing.  The  command  f2f  is  used  to  convert  the  field  file  to  the 
facet  file  to  view  the  signal  distribution  in  XCell.  The  input  parameters  and  the  response 
values  for  the  f2f  command  are  shown  in  Table  5. 
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Input  Parameters 

Response  Values 

Type  of  E-field 

Magnitude  of  Etot 

Number  of  files 

1 

Name  of  field  file 

model#real.field 

Antenna  power  level  (W) 

0.1 

Histogram  interval  (dB) 

10 

Min.  and  Max.  clip  (dBm) 

-200,  +30 

Min.  and  Max.  range  (dBm) 

-90,  +20 

Number  of  levels 

25 

Coating  code  of  the  lowest  level 

1 

Name  of  output  facet  file 

model#realout.facet 

Side  of  sample  footprint  square  (meter) 

0.5 

Shift  footprint  according  to  z-data 

Yes 

Enter  z-offset  of  footprint 

0 

Table  5.  Input  parameters  and  response  values  for  f2f  command  for  high-fidelity 

model 

Figures  23  through  32  show  the  signal  contours  for  the  high  fidelity  building 
model  with  closed  windows. 


(65  m,  -5  m)  (-10  m,  -5  m)  +  ^  dRm 


(65  m,  35  m)  L  (-10  m,  35  m).  90  dBm 


Figure  23.  WLAN  signal  contours  for  the  first  floor  for  the  parameters  in  Tables  4 

and  5 
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(65  m,  -5  m)  (-10  m,  -5  m)  +  20  dBm 


Figure  24.  WLAN  signal  contours  for  the  second  floor  for  the  parameters  in  Tables  4 

and  5 


(65  m, -5  m)  (- 1 0  m, -5  m)  +  20  dBm 


Figure  25.  WLAN  signal  contours  for  the  third  floor  for  the  parameters  in  Tables  4 

and  5 
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(65  m,  -5  m) 


(-10  m,  -5  m)+  20  dBm 


Figure  26.  WLAN  signal  contours  for  the  fourth  floor  for  the  parameters  in  Tables  4 

and  5 


(65  m,  -5  m) 


(-10  m, -5  m)  +20  dBm 


Figure  27.  WLAN  signal  contours  for  the  fifth  floor  for  the  parameters  in  Tables  4 

and  5 
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(65  m,  -5  m) 


(-10  m, -5  m)  +20  dBm 


(65  m,  35  m) 


(-10  m,  35  m)  .  90  dBm 


Figure  28.  WLAN  signal  contours  for  the  sixth  floor  for  the  parameters  in  Tables  4 

and  5 


(65  m,  -5  m) 


(-10  m,  -5  m)  +20  dBm 


Figure  29.  WLAN  signal  contours  for  the  seventh  floor  for  the  parameters  in  Tables  4 

and  5 
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(65  m,  -5  m)  (-10  m,  -5  m)  +  20  dBm 


Figure  30.  WLAN  signal  contours  for  the  eighth  floor  for  the  parameters  in  Tables  4 

and  5 


(65  m,  -5  m)  (-10  m,  -5  m)  +  20  dBm 


Figure  3 1 .  WLAN  signal  contours  for  the  ninth  floor  for  the  parameters  in  Tables  4 

and  5 
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(65  m,  -5  m) 


(-10  m,  -5  m)  +  20  dBm 


Figure  32.  WLAN  signal  contours  for  the  tenth  floor  for  the  parameters  in  Tables  4 

and  5 

As  shown  in  Figures  23  through  32,  the  plots  are  scaled  from  -90  dBm  to  +20 
dBm.  The  reason  for  this  is  to  make  the  signal  distribution  for  each  floor  more 
understandable.  Also,  a  typical  receiver  sensitivity  for  802.1  lg  is  11  Mbps  at  -90  dBm 
[21].  The  results  of  the  10  simulations  are  summarized  in  Table  6.  The  maximum  and 
minimum  signal  levels  at  each  observation  plane  inside  the  boundary  of  the  building  were 
calculated  for  the  100  mW  transmitter. 


Floor 

Minimum  Signal  Level  (dBm) 

Maximum  Signal  Level  (dBm) 

1 

-253.57 

-28.13 

2 

-  105.99 

-20.61 

3 

-  99.96 

-  13.11 

4 

-  102.86 

-5.57 

5 

-  92.02 

+  5.34 

6 

-  200.00 

+  1.37 

7 

-  200.00 

-  10.40 

8 

-  200.00 

-  18.13 

9 

-  200.00 

-  23.72 

10 

-  200.00 

-  29.33 

Table  6.  Maximum  and  minimum  signal  level  at  each  floor  with  a  WLAN 

transmitter  power  of  100  mW 
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The  minimum  and  maximum  signal  levels  for  the  fifth  floor  are  -92.02  dBm  and 
+5.34  dBm.  When  comparing  the  maximum  signal  levels,  the  reason  the  loss  at  the  sixth 
floor  is  less  than  at  the  fourth  floor  is  due  to  the  location  of  the  transmitter.  Since  the 
transmitter  is  located  at  the  ceiling  of  the  fifth  floor,  the  observation  plane  at  the  sixth 
floor  is  closer  to  the  transmitter  than  the  observation  plane  at  the  fourth  floor. 

2.  Jammer  Effects  on  the  High-Fidelity  Model 

The  second  set  of  simulations  was  run  with  the  jammer  radiating  in  the  presence 
of  the  high-fidelity  building  model.  The  properties  of  the  jammer  were  described  in  detail 
in  Section  B  of  this  chapter.  Tables  7  and  8  show  the  input  parameters  and  response 
values  for  Urbana  input  file  and  for  the  f2f  command. 


Input  Parameters 

Response  Values 

Input  file 

j  ammer  1  Otilt#.ur_input 
j  ammer20tilt#.ur_input 
j  ammer3  Otilt#.ur_input 

Facet  file 

modelfinal.facet 

Edge  file 

model0136.edge 

Antenna  file 

jammer#.antenna 

Radiation  pattern  file 

jammer#.dat 

Observation  plane  file 

obv#.txt 

Antenna  type 

10°  HPBW  directional 

20°  HPBW  directional 

30°  HPBW  directional 

Antenna  power 

1  (linear  value  for  calculation) 

Antenna  polarization 

vertical 

Antenna  frequency 

2.4  GHz 

Antenna  coordinates  (x,y,z) 

30.0,  60.0,  1.0  (meters)  and  tilted  30°  up 

Propagation  mechanism 

GO 

Edge  diffraction 

Yes  (UTD) 

Ray  angular  intervals  (degree) 

2 

Max  ray  bounces 

10 

Coating  materials 

ICOAT=  0,  1,  2,  4,  6  (as  described  in  Table  2) 

Table  7.  Input  parameters  for  the  jammer  against  the  high-fidelity  model 
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Input  Parameters 

Response  Values 

Type  of  E-field 

Magnitude  of  Etot 

Number  of  files 

1 

Name  of  field  file 

jammer  1  Otilt#.field 
jammer20tilt#.field 
j  ammer3  Otilt#.  field 

Antenna  power  level  (W) 

10  and  100 

Histogram  interval  (dB) 

10 

Min.  and  Max.  clip  (dBm) 

-200,  +30 

Min.  and  Max.  range  (dBm) 

-90,  +20 

Number  of  levels 

25 

Coating  code  of  the  lowest  level 

1 

Name  of  output  facet  file 

j  ammer  1  Otilt#.  facet 
jammer20tilt#.facet 
j  ammer3  Otilt#.  facet 

Side  of  sample  footprint  square  (meter) 

0.5 

Shift  footprint  according  to  z-data 

Yes 

Enter  z-offset  of  footprint 

0 

Table  8.  Input  parameters  and  response  values  for  f2f  command  for  the  jammer 

The  results  of  the  simulations  for  the  fifth  floor  are  summarized  in  Table  9.  When 
the  power  of  the  jammer  is  increased  by  a  factor  of  ten,  the  minimum  and  the  maximum 
signal  level  for  either  case  increases  10  dB,  as  expected. 


Jammer 

Power 

(W) 

Minimum  Signal  Level 

(dBm) 

Maximum  Signal  Level 

(dBm) 

10°  HPBW 

10 

-  88.49 

+  14.13 

100 

-  78.48 

+  24.13 

20°  HPBW 

10 

-  90.48 

+  15.20 

100 

-  80.48 

+  25.20 

30°  HPBW 

10 

-  82.09 

+  13.01 

100 

-  72.09 

+  23.01 

Table  9.  Maximum  and  minimum  signal  levels  on  the  fifth  floor  for  the  jammer 
with  three  different  HPBWs,  and  two  different  power  levels 
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3.  Signal-to-Jam  Ratio 

After  running  simulations  for  both  the  WLAN  and  jammer  independently,  the  f2fd 
command  is  then  used  to  obtain  the  differential  signal  distribution,  i.e.,  signal-to-jam  ratio 
(SJR)  for  observation  points  on  the  fifth  floor.  The  fif'd  command  can  only  be  run  by 
using  the  same  observation  plane  file  for  both  field  files.  The  input  parameters  and  the 
response  values  for  the  f2fd  command  are  shown  in  0.  Each  time  the  model5real.field  file 
is  entered  on  the  positive  (+)  side,  while  the  jammer#.field  file,  on  the  other  hand,  is 
entered  on  the  negative  (-)  side.  This  is  in  accordance  with  the  definition  of  SJR  in 
Equation  (2.33).  The  zero  point  for  the  power  difference  is  the  bum-through  condition.  If 
the  power  difference  is  greater  than  zero  (red  areas),  then  it  is  said  that  there  is  sufficient 
available  signal  for  the  receivers  on  that  particular  observation  plane  to  establish  and 
maintain  the  link. 


Input  Parameters 

Response  Values 

Type  of  E-field 

Magnitude  of  Etot 

Take  absolute  value  of  dB  difference 

No 

Number  of  field  files  on  (+)  side 

1 

Name  of  field  file  on  (+)  side 

model#real.field 

Antenna  power  level  scale  factor 

0.1 

Number  of  field  files  on  (-)  side 

1 

Name  of  field  file  on  (-)  side 

j  ammer  1  Otilt#.  field 
jammer20tilt#.field 
j  ammer3  Otilt#.  field 

Antenna  power  level  scale  factor 

10  or  100 

Histogram  interval  (dB) 

100 

Min.  and  Max.  clip  (dBm) 

-  100,  +  100 

Min.  and  Max.  range  (dBm) 

-  100,  +  100 

Number  of  levels 

2 

Coating  code  of  the  lowest  level 

1 

Name  of  output  facet  file 

m#j  10.  facet 
m#j  20.  facet 
m#j  30.  facet 

Side  of  sample  footprint  square  (meter) 

0.5 

Shift  footprint  according  to  z-data 

Yes 

Enter  z-offset  of  footprint 

0 

Table  10.  Input  parameters  and  the  response  values  for  the  f2fd  command 
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The  results  are  shown  in  Figures  33  through  38.  Each  figure  shows  the  WLAN 
signal  distribution,  the  jammer  signal  distribution,  and  the  SJR.  Only  data  for  the  fifth 
floor  is  shown,  since  this  is  the  intended  coverage  area  of  the  WLAN. 


(65  m,  -5  m) 


(-10  m,  -5  m) 


+  20  dBm 


+  100  dBm 


Figure  33.  Signal  contours  for  the  fifth  floor  (a)  100  mW  WLAN  transmitter,  (b)  10 
W  jammer  with  10°  HPBW,  (c)  Signal-to-jam  ratio  (SJR) 
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(-10  m,  -5  m) 


(65  m,  35  m) 
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X 


(-10  m,  35  m) 


Figure  34.  Signal  contours  for  the  fifth  floor  (a)  100  mW  WLAN  transmitter,  (b)  100 
W  jammer  with  10°  HPBW,  (c)  Signal- to-j am  ratio  (SJR) 
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(65  m,  -5  m) 


(-10  m,  -5  m) 
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m,  35  m) 
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-  90  dBm 


Figure  35.  Signal  contours  for  the  fifth  floor  (a)  100  mW  WLAN  transmitter,  (b)  10 
W  jammer  with  20°  HPBW,  (c)  Signal-to-jam  ratio  (SJR) 
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+  20  dBm 


(65  m,  -5  m) 


(-10  m,  -5  m) 


Figure  36.  Signal  contours  for  the  fifth  floor  (a)  100  mW  WLAN  transmitter,  (b)  100 
W  jammer  with  20°  HPBW,  (c)  Signal-to-j am  ratio  (SJR) 
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(65  m,  35  m) 
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-  90  dBm 


+  100  dBm 


Figure  37.  Signal  contours  for  the  fifth  floor  (a)  100  mW  WLAN  transmitter,  (b)  10 
W  jammer  with  30°  HPBW,  (c)  Signal-to-jam  ratio  (SJR) 
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-  90  dBm 


Figure  38.  Signal  contours  for  the  fifth  floor  (a)  100  mW  WLAN  transmitter,  (b)  100 
W  jammer  with  30°  HPBW,  (c)  Signal-to-jam  ratio  (SJR) 
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As  shown  in  Figures  33  through  38,  for  all  cases,  the  jammer  mostly  affects  the 
rooms  closest  to  the  jammer.  In  part  (b)  of  each  figure,  due  to  the  shadow  regions  at  the 
comer  of  the  building,  the  jammer  effect  decreases  rapidly.  The  rooms  on  the  far  side  of 
the  building  are  less  affected  when  the  jammer  is  10  W.  When  the  power  of  the  jammer  is 
increased  to  100  W,  the  available  WLAN  signal  area  then  obviously  decreases  on  both 
sides  of  the  building.  In  order  to  quantify  the  results,  the  number  of  blue  and  red  points 
are  calculated  and  written  as  a  percentage.  The  percentage  value  of  the  available  signal  in 
part  (c)  of  each  figure  is  also  tabulated.  Table  1 1  shows  the  results  of  the  simulations. 


Floor 

Jammer 

Power 

Difference  (dB) 

#  of  Points 

Available 
Signal  in 
the  Area 

HPBW 

Power  (W) 

Min 

Max 

(-)  Side 

(+)  Side 

5 

10° 

10 

-81.53 

64.79 

1026 

2705 

72.50% 

100 

-91.53 

54.79 

1917 

1814 

48.62% 

20° 

10 

-82.55 

54.93 

1250 

2481 

66.50% 

100 

-92.55 

44.93 

2196 

1535 

41.14% 

30° 

10 

-80.10 

54.41 

1532 

2199 

58.94% 

100 

-90.10 

44.41 

2598 

1133 

30.37% 

Table  1 1 .  The  results  of  the  f2fd  command  for  the  fifth  floor 


The  percentage  values  in  Table  1 1  show  that  as  either  the  HPBW  or  the  power  of 
the  jammer  increases,  the  available  WLAN  signal  in  the  area  decreases.  The  results  for 
the  variations  in  HPBW  and  power  are  illustrated  in  Figure  39.  The  Matlab  code  to  plot 
Figure  39  appears  in  the  Appendix. 
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Figure  39.  The  effects  of  different  HPBWs  and  power  levels  for  the  fifth  floor 

E.  SUMMARY 

This  chapter  began  with  the  description  of  the  WLAN  transmitting  antenna,  the 
jammer  and  the  observation  plane,  and  then  proceeded  to  study  the  simulations.  The 
simulations  were  conducted  for  three  scenarios.  First,  the  signal  distribution  of  a  100  mW 
WLAN  transmitting  antenna  was  examined  for  all  floors.  Then,  the  jammer  effects, 
including  three  different  HPBWs  with  two  different  power  levels,  were  examined  only 
for  the  fifth  floor.  Finally,  the  signal-to  jam  ratio  (SJR)  was  computed  for  the  fifth  floor. 
The  results  for  other  floors  are  tabulated  and  shown  in  the  Appendix.  The  next  chapter 
includes  a  summary,  conclusions,  and  recommended  future  work  relating  to  this  research. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  SUMMARY 

The  primary  objective  of  this  thesis  was  to  investigate  the  performance  of  wireless 
local  area  networks  (WLANs)  in  an  indoor  environment  with  the  presence  of  noise 
jamming  and  calculate  the  signal-to-jam  ratio  (SJR)  in  a  high-rise  building.  The  high-rise 
building  model  was  developed  with  the  Rhino  software  and  the  Urbana  software  was 
used  to  calculate  the  signal  levels  for  the  simulations.  Specifically,  the  frequency  of  2.4 
GHz  was  studied.  In  order  to  make  the  entire  building  physically  secure,  a  low  power 
(100  mW)  WLAN  transmitter  was  located  at  the  ceiling  of  the  fifth  floor.  The  jammer 
was  located  30  meters  away  from  the  building  in  the  parking  area.  The  real  conditions 
were  considered  for  the  simulations  and  SJR  was  only  calculated  for  the  area  inside  of  the 
high-fidelity  building  model.  The  results  indicate  that  increasing  either  the  HPBW  or  the 
power  of  the  jammer  reduces  the  available  WLAN  signal. 

B.  CONCLUSIONS 

The  Urbana  software  is  very  useful  tool  for  modeling  wireless  signal  distribution 
for  indoor  environments.  Signal  coverage  and  losses  can  be  analyzed  easily,  and  the 
effectiveness  of  the  WLAN  can  be  optimized  by  running  Urbana  simulations.  If  the 
signal  coverage  exceeds  the  area  where  needed,  i.e.,  WLAN  signals  propagate  from  the 
inside  of  the  building  to  the  outside,  some  receivers  with  high  sensitivities  might  intercept 
the  leaking  WLAN  signals.  With  the  help  of  simulations,  the  security  posture  of  the 
WLAN  can  be  improved  as  follows, 

•  decreasing  the  power  of  the  WLAN  transmitter, 

•  increasing  the  sensitivity  of  the  access  point  and  the  receivers, 

•  locating  the  transmitting  antenna  in  the  most  interior  part  of  the  building, 
and 

•  using  directional  antennas  to  limit  the  signal  distribution. 

If  the  WLAN  signal  is  not  interceptable  due  to  limited  coverage  and  physical 
security,  then  the  alternative  is  to  jam  the  signal.  The  effectiveness  of  the  jammer 
depends  on, 
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•  the  distance  to  the  WLAN  coverage  area, 

•  the  power  of  the  jammer  transmitter, 

•  the  radiation  pattern  of  the  jammer  antenna, 

•  the  material  properties  of  the  building,  and 

•  the  openings  of  the  building. 

C.  FUTURE  WORK 

The  present  research  was  only  concerned  with  calculating  S  JR  for  a  single  WLAN 
transmitter  versus  a  single  jammer.  Future  work  might  entail, 

•  analyzing  the  SJR  for  multiple  WLAN  transmitters  versus  multiple 
jammers, 

•  comparing  the  effects  of  different  placements  of  the  transmitters,  such  as 
placing  the  jammer  inside  the  building. 

•  analyzing  the  effects  of  different  polarizations  of  the  WLAN  transmitter 
antenna  and  the  jammer  antenna, 

•  analyzing  the  SJR  by  using  different  types  of  antennas, 

•  analyzing  the  effects  of  different  material  properties,  such  as  using  PEC 
materials  for  the  outer  walls,  and 

•  calculating  the  jam-to-signal  ratio  (JSR). 
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APPENDIX 


The  appendix  includes  a  sample  Urbana  input  file  with  and  antenna  file,  Matlab 
codes  for  the  jammer  antenna,  the  observation  points  file  and  the  SJR.  Finally,  the 
tabulated  values  and  the  plots  are  also  included  for  the  other  floors. 

A.  SAMPLE  URBANA  INPUT  FILE  AND  ANTENNA  FILE 

-  input  Urbana  v  2 . 5 - 

# 

jj:  ******************************** 

#  A - scatterer  file, length  &  freq 

jj:  ******************************** 

#  -  name  of  scatterer  file  in  ACAD  format  (e.g.  wall. facet) 

modelf inal . facet 

#  -  length  unit : l=inch,  2=cm,  3=meter,  4=mm,  5=mil 

3 

#  -  uniform  freq  (GHz) :  start  freq,  end  ,  nstep 

#  (nstep=0  means:  just  do  first  freq.  CAUTION:  antenna  patterns  are 

#  assumed  to  be  indep.  of  freq  and  is  calculated  at  end  freq) 

2.4  2.4  0 

# 

jj:  ********************************* 

#  B -  Antenna  Description  and  List 

jj:  ********************************* 

# 

#  - Enter  method  of  describing  antennas. 

#  (1  =  here,  2  =  file) : 

2 

#  - If  described  in  file,  enter  file  name: 

jammertilt . antenna 

#  - If  described  here,  fill  in  sections  Bl,  B2,  B3 . 

#  If  described  in  file,  use  dummy  data  in  sections  Bl,  B2,  B3 

#  (specify  one  dummy  antenna  type,  dummy  antenna  origin, 

#  and  one  dummy  item  in  antenna  list) . 

# 

jj:  ************************ 

#  Bl :  Define  Antenna  Types 
ft  ************************ 

# 

#  Two  lines  for  each  type. 

#  Linel:  type  ID,  ant  code 

#  Line2 :  parameters 

# 

#  Type  ID  must  start  from  1  and  increment  by  1  thereafter 

# 

#  Ant  Code  meaning  parameters 

# 

#  1  pattern  file  filename (ascii ) 

#  2  dipole  length (real) 

# 

#  Antenna  Types  list: 

# 

#  Enter  number  of  antenna  types : 

1 

#  Type  #1 
1  1 

0.0625 

# 

§  *********************************************** 

#  B2 :  Enter  origin  of  antenna  coord  in  main  coord 
§  *********************************************** 

# 

0.00.0.  .0 

# 
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§  *********************** 

#  B3 :  Create  Antenna  List 
§  *********************** 

# 

#  Three  lines  for  each  antenna. 

#  Linel:  Type  ID,  location  (x,y,z),  power  (watts),  phase (deg) 

#  Line2 :  Local  x-axis  in  main  coord. 

#  Line3:  Local  z-axis  in  main  coord. 

# 

#  Enter  number  of  antennas : 

1 

# 

#  Antenna  #1 

1  30.0  60.0  1.0  1  0 

-1.  0.5  0.8660 

0.  -0.8660  0.5 

# 

jj:  ********************** 

#  C - Observation  points 

jj:  ********************** 

#  -  Observation  points  defined  with  respect  to  main  coord,  system  7. 

#  Enter  method  of  specifying  list  of  points. 

#  (1  =  here,  2  =  file) : 

2 

#  -  If  points  are  listed  here,  enter  number  of  points  (kobtot) : 

1 

#  -  If  listed  here  (1  above).  List  xyz  of  points  in  main  coord  7 

#  (one  point  at  a  line) .  If  2  above,  include  one  dummy  line. 

1.  2.  -11.00 

#  -  If  points  listed  in  file  (2  above),  enter  name  of  file. 

obv5 . txt 

#  -  Include  direct  Tx  to  observer  contribution. 

#  If  you  turn  on  the  direct  contribution  from  the  transmitter  to  the 

#  observation  point,  computed  result  will  be  the  total  field,  which  is 

#  the  incident  +  scattered  field.  For  propagation  analysis,  this  is 

#  the  preferred  setting.  Otherwise,  the  result  only  includes  the 

#  scattered  field. 

# 

#  Include  direct  contribution  from  transmitter  to  observation  point  (rx) 

#  (1  =  yes,  0,2  -  no) : 

1 

#  -  Compute  received  power  into  Rx  antenna. 

#  Urbana  always  computes  field  levels  at  the  observation  point. 

#  If  you  specify  an  Rx  antenna,  Urbana  will  also  compute  the  received 

#  power  and  record  the  results  in  the  (runname) .couple  file. 

#  This  causes  a  moderate  but  slow-down  when  using  the  SBR  method  (below) . 

# 

#  Include  Rx  antenna  (1  =  yes,  0,2  =  no) : 

0 

#  -  Rx  antenna  specification 

#  Remaining  entries  in  Section  C  can  be  ignored  if  not  including 

#  an  Rx  antenna. 

#  Enter  antenna  type  (1  =  pattern  file,  2  =  dipole) : 

2 

#  Each  antenna  type  requires  additional  parameters. 

#  List  of  expected  parameters  follows.  Choose  one. 

# 

#  Type  Description  Expected  Parameter (s) 

#  1  Pattern  File  File  Name  (e.g.,  beam.antpat) 

#  2  Dipole  Length  (in  prevailing  unit) 

# 

#  Enter  parameter (s)  on  next  line: 

2.5 

#  -  Rx  antenna  orientation 

#  Enter  local  x-axis  of  Rx  in  global  coordinates 

1.  0.  0. 

#  Enter  local  z-axis  of  Rx  in  global  coordinates 

0.  0.  1. 

# 

jj:  ***************************** 

#  D - Theoretical  consideration 

jj:  ***************************** 


56 


#  -  Choose  method  of  computation 

#  0  =  compute  fields  in  the  ABSENCE  of  the  scatterer 

#  1  =  compute  fields  by  SBR 

#  2  =  compute  fields  by  GO 

2 

#  -  If  SBR,  select  a  PO  integration  scheme  at  bounce  points 

#  1  =  do  integration  at  first  &  last  bounce  points  only 

#  2  =  do  so  at  all  bounce  points  (GTD  formulation) 

#  3  =  do  so  at  all  bounce  points  (PTD  formulation) 

2 

#  -  Edge  diffraction 

#  SBR  can  be  enhanced  with  PTD  edge  diffraction. 

#  GO  can  be  enhanced  with  GTD  edge  diffraction. 

#  Add  edge  diffraction  (0,2=no,  1=ILDC  (SBR  or  GO),  3=UTD  (GO  only) 

3 

#  -  If  edge  diffraction  switched  on,  enter  name  of  edge  file 

#  (e.g.,  wall. edge  or  dummy  if  edge  not  included) . 
model0136 . edge 

#  -  Choose  method  of  ray  launch 

#  1  =  by  (baby)  facet,  achieving  a  uniform  first  bounce  surface  density 

#  2  =  uniform  angular  distribution  (burst  launch) 

#  (If  computation  by  GO,  must  select  2  =  burst  launch) 

2 

#  -  If  ray  launch  by  (baby)  facet  (1  above),  enter  ray  density: 

#  #  rays/wavelength  (normally  5-10) 

5 

#  -  If  burst  ray  launch  (2  above) ,  enter  angular  interval  (deg) . 

#  (Typically  0.25  -  2.0  deg) 

2. 

#  -  max  permissible  ray  bounces  (normally  5-10) 

10 

#  -  max-voxdepth  =  max  depth  of  BSP  tree  (normally  20) 

#  max-voxl  =  max  facets  in  each  voxel (normally  10) 

#  (Larger  voxdepth  &  smaller  voxl  lead  to  faster  ray  tracing 

#  but  more  computer  memory) 

20,10 

#  -  ICOAT  for  absorbing  facets 

888 

#  -  IQMATRIX  for  divergence  factor 

#  1  =  calculated  by  Q-matrix 

#  2  =  ignored  except  for  the  spherical  wave  spread 

2 

#  -  IF  using  Q-matrix,  name  target  curvature  file (e.g.  wall.curv) 

dummy . curv 

#  -  IPEC=1  if  all  pec,  =2  if  coating  present 

2 

#  -  For  PEC  scatterer,  give  the  magnitude  of  reflection  coeff 

#  (use  1.0  for  ideal  PEC,  use  less  for  rough  PEC — fudging) 

1 . 

#  -  IF  PEC,  the  rest  coating  info  is  dummmy 

#  -  material  reflection  is  done  through  a  look-up  table 

#  specify  the  freq  interval  in  GHz  for  the  table  e.g.  0.25 

#  (dummy  if  input  freq  less  than  51) 

0.2 

******************** 

E - coating  material 

******************** 
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# 


—  number  of  materials 

(NOT  including  pec,  which  is  identified  by  ICOAT=0) 
(NOT  including  absorbing  facets:  ICOAT=28  or  888) 
(If  3  material,  urbana  reads  only  ICOAT=l-3) 


< - NCOTOT 

for  each  material. 


identify  its  boundary  type: 


iboundary  =  1 
2 

3 

4 

5 

6 
7 

for  each  material, 
A  ICOAT=l  AA^ 


if  impedance  boundary 

if  layered  slabs  with  PEC  backing 

if  penetrable  layered  slabs  with  air  backing 

if  penetrable  FSS  Tx/Refl  table  supplied  by  user 

if  same  as  2  except  using  freq-dep  ramlib.d 

if  antenna  refl  table  supplied  by  user 

if  layers  over  semi-infinite  ground 

given  info  by  following  the  templates 
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-  iboundary 

3 

-  number  of  layers  over  air  backing 

(1st  layer  is  farthst  fr  incid  field  and  innermost) 

1 

-  thick, epsilon (c) , mu (c)  , resistivity (ohm) 

.300  (10. 1,-0. 5)  (1.0, -0.0)  le30 

AAA  ICOAT=2  aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa. 

-  iboundary 

3 

-  number  of  layers  over  air  backing 

(1st  layer  is  farthst  fr  incid  field  and  innermost) 

1 

-  thick, epsilon (c) , mu (c) , resistivity (ohm) 

.100  (10. 1,-0. 5)  (1.0, -0.0)  le30 

AAA  ICOAT=4  ^aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa. 

-  iboundary 

3 

-  number  of  layers  over  air  backing 

(1st  layer  is  farthst  fr  incid  field  and  innermost) 

1 

-  thick, epsilon (c) , mu (c) , resistivity (ohm) 

.0250  (6.0, -0.0)  (1.0, -0.0)  le30 

AAA  ICOAT=6  aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa. 

-  iboundary 

7 

-  number  of  layers  over  air  backing 

(1st  layer  is  farthst  fr  incid  field  and  innermost) 

1 

-  thick, epsilon (c) , mu (c) , resistivity (ohm) 

.1000  (3.0, -0.0)  (1.0, -0.0)  le30 

- epsilon (c) , mu (c)  of  semi-infinite  ground 

(3.00,-0.00)  (1.00,-0.00) 


(End  of  regular  input  file.  Leave  a  few  blank  lines) 


’OPTIONAL  ADVANCE  FEATURES’  (Do  not  change  letters  in  quotations) 

#  The  line  above  must  be  placed  at  the  end  of  the  regular  urbana 

#  input.  Advance  features  are  designed  for  special  applications  or 

#  for  testing  codes.  They  are  not  needed  by  general  usages. 

#  - 

#  ADVANCE 1 :  ADD  GTD-TYPE  BLOCKAGE  CHECK 

#  - 

#  In  regular  urbana  computation,  blockage  check  is  mostly  done  by 

#  PTD  principle.  For  interior  scattering  in  a  confined  region,  use  of 

#  GTD  principle  may  be  more  appropriate. 

#  Option  to  use  GTD  principle:  l=yes,  2=no  (regular  case) 

2 

#  - 

#  ADVANCE2:  SIMPLE  TERRAIN  BLOCKAGE  MODEL 

#  - 

#  For  GO  method,  terrain  generates  100%  blockage,  and  blocked  rays  leave 

#  no  energy  behind  a  hill.  With  this  feature,  LOS  rays  and  UTD  edge 

#  diffraction  rays  can  pass  through  terrain,  with  some  attenuation. 

#  Attenuation  is  measured  in  dB  per  hill.  Each  hill  is  identified 

#  by  two  passages  through  two  terrain  facets. 

#  Can  only  be  used  with  GO  method  (and  UTD  edge  option) . 

#  Use  simple  terrain  model:  1  =  yes,  2  =  no  (regular  case) 

2 

#  Enter  coating  code  range  of  terrain  facets  (e.g.,  1,  2) : 

11 

#  Enter  amount  of  attenuation  per  hill  (dB,  >  0) : 

5. 

# 

#  ADVANCE 3 :  APPROXIMATE  DOUBLE  DIFFRACTION  MODEL 

#  - 

#  For  GO  +  UTD  method,  only  single  diffraction  is  considered. 

#  With  this  feature,  double  diffraction  is  approximated  by  identifying 
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#  surfaces  which  block  the  single  diffraction,  such  as  building  walls. 

#  If  one  or  two  facets  block  the  path  from  the  single  diffraction  point 

#  to  the  transmitter,  the  diffraction  is  still  included,  but  with  attenuation. 

#  Works  best  if  "diffracting  facets",  marked  by  their  coating  code,  are 

#  always  associated  with  enclosed  structures  with  well  defined  edges. 

#  Use  double  diffraction  model:  1  =  yes,  2  =  no  (regular  case) 

2 

#  Encounter  coating  code  range  of  diffracting  facets  (e.g.,  5,  10)  : 

2  2 

#  Enter  amount  of  attenuation  for  second  diffraction  (dB,  >  0); 

10. 

#  - 

#  ADVANCE 4 :  ACCELERATION 

#  - 

#  For  large  scenes,  run  time  grows  both  with  the  number  of  field 

#  observation  points  and  the  number  of  edges.  Normally,  all  combinations 

#  of  lit  edges  and  observation  points  are  considered.  This  feature 

#  accelerates  the  processing  by  limiting  the  scope  of  considered  edge 

#  interactions  to  region  around  the  LOS  path  from  the  transmitter 

#  to  the  observation  point.  For  example,  to  run  a  5  km  by  5  km  scene, 

#  one  may  choose  a  250  m  interaction  radius.  For  each  observation 

#  point,  edges  are  ignored  that  lie  outside  an  ellipse  whose  foci  are  the 

#  Tx  and  the  observation  point  and  whose  major  axis  is  the  LOS  distance 

#  plus  500  m  (radius  x  2) . 

#  This  feature  can  also  be  used  to  automatically  filter  edge  files 

#  whose  domain  far  exceeds  the  domain  of  observation  points. 

#  Only  use  this  feature  for  terrestrial  simulations  where  the  scene 

#  is  nominally  parallel  to  the  x-y  plane. 

# 

#  Use  large  scene  acceleration:  1  =  yes,  2  =  no  (regular  case) 

2 

#  Enter  radius  of  interaction 
250. 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


************************ 

Bl :  Define  Antenna  Types 
************************ 

Two  lines  for  each  type. 

Linel:  type  ID,  ant  code 
Line2 :  parameters 

Type  ID  must  start  from  1  and  increment  by  1  thereafter 
Ant  Code  meaning  parameters 

1  pattern  file  filename (ascii ) 

2  dipole  length (real) 

Antenna  Types  list: 

Enter  number  of  antenna  types: 

Type  #1 

1 

ammerl . dat 

*********************************************** 

B2 :  Enter  origin  of  antenna  coord  in  main  coord 
*********************************************** 

.00.0.  .0 
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#  *********************** 

#  B3 :  Create  Antenna  List 
§  *********************** 

# 

#  Three  lines  for  each  antenna. 

#  Linel:  Type  ID,  location  (x,y,z),  power  (watts),  phase (deg) 

#  Line2 :  Local  x-axis  in  main  coord. 

#  Line3:  Local  z-axis  in  main  coord. 

# 

#  Enter  number  of  antennas : 

1 

# 

#  Antenna  #1 

1  30.0  60.0  1.0  1  0 

-1  0.5  0.8660 
0  -0.8660  0.5 

# 


B.  MATLAB  CODES 

1.  Jammer  Antenna  Radiation  Pattern  File 


%  Antenna  Pattern  File  Generation 

close  all; 

clc; 

clear; 

lamda  =  (3*10A8) / (2400*10A6) ; 

HPBW  =10;  %  in  degrees 

a  =  58.4  *  lamda  /  (2*HPBW)  %  Uniform  Circular  Aperture  Size  (radius) 
beta  =  2*pi/lamda; 

E0  =  1; 

c  =  j *beta*E0*pi*aA2/ (2*pi) ;  %Constant  term  in  E  field 

Directivity  =  ( 4 *pi*pi*aA2 ) /lamdaA2 

Direct ivity_dB  =  10*logl0 (Directivity) 

stepl  =  90;  %Vertical  Steps 

step2  =  180;  %Horizontal  Steps 

d_theta  =  180/ (stepl); 

d_phi  =  360/(step2); 

k  =  0; 

ip  =  0; 

iq  =  0; 

for  phi  =  0:d_phi:360 

for  theta  =  0 : d_theta : 180  %plots  out  the  zero  point 
k  =  k+1; 

arg=beta*a*sin (theta*pi/180 ) ; 

f  =  1; 

if  abs (arg) >=le-5 

f  =  2 ^bessel j ( 1 , arg) /arg; 

end 

e_theta  =  cos  (phi*pi/180 ) ^c^f ; 

e_phi  =  -sin (phi^pi/180 ) *cos (theta^pi/180 ) ^c^f ; 

if  theta  >  90 

e_theta  =  0; 

e_phi  =  0; 

end 

if  phi  ==  0 
ip  =  ip  +  1; 

etO(ip)  =  20*logl0 (abs (e_theta) ) ; 
epO(ip)  =  20*logl0 (abs (e_phi) ) ; 
eth(ip)  =  theta; 
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end 

if  phi  ==  90 
iq  =  iq  +  1; 

et90(iq)  =  20*logl0 (abs (e_theta) ) ; 
ep90 (iq)  =  20*logl0 (abs (e_phi) ) ; 
eth90(iq)  =  theta; 
end 

A(k,l:4)  =  [real (e _ theta) ,  imag (e_theta) ,  real(e_phi),  imag (e_phi ) ] ; 

end 

end 

save  antennapatl  A  -ASCII; 
figure ( 1 ) 

plot (eth, etO, eth, epO) 

title  (  1  Directional  Antenna  Pattern  (HPBW=10  degrees)  ') 
xlabel ( 'Angle  [degree] ' ) 
ylabel (' Relative  Power  [dB] ' ) 
grid 

figure  (2 ) 

plot (eth90, et90, eth90,  ep90) 
xlabel ( 'Angle  [degree] ' ) 
ylabel (' Relative  Pattern,  dB ' ) 
grid 


2.  Observation  Points  File 


%  Observation  Points 

clc; 

clear; 

step  =  .5;  %footprint 

i  =  i; 

z  =  15; 

for  x  =  -10:step:65 
for  y  =  -5:step:35 
M (i,  :  )  =  [x,  y, z]  ; 
i=i+l ; 

end 

end 


Lines=i-1 


save  observlstrip . obv  M  -ASCII; 


3.  Plotting  the  Effects  of  Different  HPBWs  and  Power  Levels 


clear;  clc; 
width= . 4 ; 

power={'10  W',  '100  W'}; 

x= [72.50,  48.62]; 
y= [66.50,  41.14]; 
z= [ 58 . 94 ,  30.37]; 

subplot  (1,3,1) 

bar  (x,  width,  'r',  'EdgeColor',  'r') 
text (0.6, max (x) +5, [ ' \bf ' , num2str (x ( 1 ) ) , '%']) 
text ( 1 . 6, min (x) +5, [ ' \bf ' , num2str (x (2 ) ) , '%']) 
axis ( [ 0  3  0  100 ]) ; 
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set (gca, ' xticklabel ' , power ,  'FontSize' ,  8); 

xlabel  (  ' \bf  10Ao  HPBW ' ,  ' FontSize ' ,  10); 

ylabel  (  ' \bf  Percentage ' , 'FontSize’,  10); 

subplot  (1,3,2) 

bar (y,  width,  'r',  'EdgeColor',  'r') 
text ( 0 . 6, max (y) +5, [ ' \bf ' , num2str (y(l)), '  %  '  ]  ) 
text (1.6, min (y) +5, [ ' \bf ' , num2str (y (2  )  )  , ’%']) 
set (gca,  ’ xticklabel ’ , power,  ’ FontSize ' ,  8 )  ; 

xlabel  (  ' \bf  20Ao  HPBW FontSize ’ ,  10); 

axis ([030  100] ) ; 

title ( { 1 \bf  Available  Signal  Coverage  at  5A{th} 

’ } ,  ’ FontSize ' ,  10 ) ; 

subplot  (1,3,3) 

bar  (z,  width,  ’r’,  'EdgeColor’,  'r') 
text (0.6, max ( z ) +5, [ ' \bf ' , num2str (z(l)), '%']) 
text ( 1 . 6, min ( z ) +5, [ ' \bf ' , num2str (z(2)), '%']) 
set (gca,  ' xticklabel ' , power,  ' FontSize ' ,  8 )  ; 

xlabel  (  ' \bf  30Ao  HPBW FontSize ' ,  10); 
axis  (  [ 0  3  0  100 ])  ; 


Floor ' , ' 


62 


c. 


TABULATED  RESULTS  AND  PLOTS  FOR  THE  OTHER  FLOORS 


Floor 

Jammer 

Power  Difference 
(dB) 

#  of  Points 

Available 
Signal  in 
the  Area 

HPBW 

Power 

(W) 

Min 

Max 

(-)  Side 

(+)  Side 

1 

10° 

10 

-181.79 

13.08 

6718 

53 

0.79% 

100 

-191.79 

3.08 

6769 

2 

0.03% 

20° 

10 

-191.51 

8.57 

6751 

20 

0.30% 

100 

-201.51 

-1.42 

6771 

0 

0.00% 

30° 

10 

-187.57 

9.77 

6745 

26 

0.38% 

100 

-197.57 

-0.23 

6771 

0 

0.00% 

Table  12.  The  results  of  the  f2fd  command  for  the  first  floor 


Available  Signal  Coverage  at  1st  Floor 


10°  HPBW  20°  HPBW  30°  HPBW 

Figure  40.  The  effects  of  different  HPBWs  and  power  levels  for  the  first  floor 
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Floor 

Jammer 

Power  Difference 
(dB) 

#  of  Points 

Available 
Signal  in 
the  Area 

HPBW 

Power 

(W) 

Min 

Max 

(-)  Side 

(+)  Side 

2 

10° 

10 

-73.36 

31.46 

2672 

1059 

28.38% 

100 

-83.35 

21.47 

3461 

270 

7.24% 

20° 

10 

-75.69 

32.47 

2634 

1097 

29.40% 

100 

-85.69 

22.47 

3478 

253 

6.78% 

30° 

10 

-83.60 

31.72 

3125 

606 

16.24% 

100 

-93.60 

21.72 

3656 

75 

2.01% 

Table  13.  The  results  of  the  f2fd  command  for  the  second  floor 
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Available  Signal  Coverage  at  2nd  Floor 
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29.4% 
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10W  100W 
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30°  HPBW 


Figure  41.  The  effects  of  different  HPBWs  and  power  levels  for  the  second  floor 
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Floor 

Jammer 

Power  Difference 
(dB) 

#  of  Points 

Available 
Signal  in 
the  Area 

HPBW 

Power 

(W) 

Min 

Max 

(-)  Side 

(+)  Side 

3 

10° 

10 

-77.56 

49.34 

1737 

1994 

53.44% 

100 

-87.56 

39.34 

2807 

924 

24.77% 

20° 

10 

-79.63 

59.08 

1825 

1906 

51.09% 

100 

-89.63 

49.08 

2975 

756 

20.26% 

30° 

10 

-83.07 

37.95 

2697 

1034 

27.71% 

100 

-93.07 

27.95 

3464 

267 

7.16% 

Table  14.  The  results  of  the  f2fd  command  for  the  third  floor 

Available  Signal  Coverage  at  3rd  Floor 


10°  HPBW  20°  HPBW  30°  HPBW 


Figure  42.  The  effects  of  different  HPBWs  and  power  levels  for  the  third  floor 
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Floor 

Jammer 

Power  Difference 
(dB) 

#  of  Points 

Available 
Signal  in 
the  Area 

HPBW 

Power 

(W) 

Min 

Max 

(-)  Side 

(+)  Side 

4 

10° 

10 

-80.01 

127.05 

1563 

2168 

58.11% 

100 

-90.01 

127.05 

2599 

1132 

30.34% 

20° 

10 

-83.96 

127.05 

1768 

1963 

52.61% 

100 

-93.96 

127.05 

2896 

835 

22.38% 

30° 

10 

-85.00 

127.05 

2339 

1392 

37.31% 

100 

-95.00 

127.05 

3300 

431 

11.55% 

Table  15.  The  results  of  the  f2fd  command  for  the  fourth  floor 

Available  Signal  Coverage  at  4th  Floor 


10°  HPBW  20°  HPBW  30°  HPBW 


Figure  43.  The  effects  of  different  HPBWs  and  power  levels  for  the  fourth  floor 
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Floor 

Jammer 

Power  Difference 
(dB) 

#  of  Points 

Available 
Signal  in 
the  Area 

HPBW 

Power 

(W) 

Min 

Max 

(-)  Side 

(+)  Side 

6 

10° 

10 

-123.97 

128.52 

1545 

2186 

58.59% 

100 

-133.97 

128.52 

2543 

1188 

31.84% 

20° 

10 

-121.61 

128.52 

1833 

1898 

50.87% 

100 

-131.61 

128.52 

2903 

828 

22.19% 

30° 

10 

-127.20 

128.52 

2122 

1609 

43.12% 

100 

-137.20 

128.52 

3122 

609 

16.32% 

Table  16.  The  results  of  the  f2fd  command  for  the  sixth  floor 
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Figure  44.  The  effects  of  different  HPBWs  and  power  levels  for  the  sixth  floor 
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Floor 

Jammer 

Power  Difference 
(dB) 

#  of  Points 

Available 
Signal  in 
the  Area 

HPBW 

Power 

(W) 

Min 

Max 

(-)  Side 

(+)  Side 

7 

10° 

10 

-113.61 

127.84 

2094 

1637 

43.88% 

100 

-123.61 

127.84 

3137 

594 

15.92% 

20° 

10 

-117.04 

127.84 

2497 

1234 

33.07% 

100 

-127.4 

127.84 

3416 

315 

8.44% 

30° 

10 

-119.49 

127.84 

2640 

1091 

29.24% 

100 

-129.49 

127.84 

3562 

169 

4.53% 

Table  17.  The  results  of  the  f2fd  command  for  the  seventh  floor 


Available  Signal  Coverage  at  7th  Floor 


10°  HPBW  20°  HPBW  30°  HPBW 

Figure  45.  The  effects  of  different  HPBWs  and  power  levels  for  the  seventh  floor 
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Floor 

Jammer 

Power  Difference 
(dB) 

#  of  Points 

Available 
Signal  in 
the  Area 

HPBW 

Power 

(W) 

Min 

Max 

(-)  Side 

(+)  Side 

8 

10° 

10 

-120.16 

139.07 

2799 

932 

24.98% 

100 

-130.16 

139.07 

3471 

260 

6.97% 

20° 

10 

-122.50 

139.07 

3160 

571 

15.30% 

100 

-132.50 

139.07 

3590 

141 

3.78% 

30° 

10 

-130.26 

139.07 

3317 

414 

11.10% 

100 

-140.26 

139.07 

3700 

31 

0.83% 

Table  18.  The  results  of  the  f2fd  command  for  the  eighth  floor 
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Figure  46.  The  effects  of  different  HPBWs  and  power  levels  for  the  eighth  floor 
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Table  19.  The  results  of  the  f2fd  command  for  the  ninth  floor 
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Figure  47.  The  effects  of  different  FIPBWs  and  power  levels  for  the  ninth  floor 
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Table  20.  The  results  of  the  fif'd  command  for  the  tenth  floor 
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Figure  48.  The  effects  of  different  HPBWs  and  power  levels  for  the  tenth  floor 
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